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ABSTRACT 

We construct time-dependent models of accretion discs around black holes and neutron 
stars. We investigate the effect of evaporating the disc inner regions during quiescence 
on the predictions of the Disc Instability Model (DIM) for these systems. We do not 
include irradiation of the disc in the models. 

Removing the inner, most unstable parts of the accretion disc increases the 
predicted recurrence times. However, DIMs with values of the viscosity parameter 
tthot ~ 0.1 and Q!coid ^ 0.02 (typical of applications of the DIM to standard dwarf 
nova outbursts) fail to reproduce the long recurrence times of SXTs (unless we resort 
to fine-tuning of the parameters) independent of the evaporation strength. We show 
that models with evaporation and a smaller value of acoid 0.005) do reproduce 
the long recurrence times and the accretion rates at the level of the Eddington rate 
observed in outburst. The large difference between the values of Ofhot and acoid, if con- 
firmed when disc irradiation is included, suggests that several viscosity mechanisms 
operate in these accretion discs. 

For some choices of parameters our models predict reflares during the decline 
from outburst. They are a physical property of the model and result from a heating 
front forming in the wake of an initial cooling front and subsequent, multiple front 
reflections. The reflares disappear in low-a models where front reflection can not 
occur. 

Key words: X-ray: stars - binaries: close - accretion, accretion discs - black hole 
physics- novae, cataclysmic variables - instabilities - turbulence - MHD 



1 INTRODUCTION 

Low Mass Black Hole Binaries (LMBHBs) and Low Mass 
Neutron Stars Binaries (LMNSBs) are close binary systems 
in which a Roche-lobe filling, low mass (main sequence or 
sub-giant) secondary star transfers mass to a compact pri- 
mary. All known LMBHBs and several LMNSBs are tran- 
sients (see e.g. Tanaka & Lewin 1995, van Paradijs & Mc- 
Clintock 1995, White, Nagase & Parmar 1995); these tran- 
sient systems are known as soft X-ray transients (SXTs), or 
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X-ray novae. Their outbursts last from several weeks to sev- 
eral months during which X-ray luminosities can approach 
the Eddington limit. The outbursts are separated by long 
periods of quiescence lasting from one to tens of years (see 
e.g. Chen, Shrader & Livio 1997; hereafter CSL). 

SXTs are in some respects similar to dwarf novae (DN), 
which exhibit outbursts of 4 - 6 magnitudes in optical. These 
outbursts typically last days and are separated by periods of 
quiescence of a few weeks. DN are a subclass of Cataclysmic 
Variables (CVs; see Warner 1995 for a review), close binary 
systems similar to LMBHBs and LMNSBs, except for the 
compact primary which is a white dwarf. 

Is is widely believed that dwarf nova outbursts are due 
to a thermal- viscous instability in the accretion disc of these 
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systems. For effective temperatures between about 5000 K 
and 8000 K, accretion discs in which viscosity is described by 
the "a— prescription" (Shakura & Sunyaev 1973) are ther- 
mally and viscously unstable because of large changes in 
the opacity when hydrogen recombines. The disc instabil- 
ity model (DIM; see Cannizzo 1993b for a review; Hameury 
et al. 1998 and Menou, Hameury & Stehle 1999a for the 
most recent version of the model) accounts for a number of 
general properties of dwarf novae. Similarities between the 
properties of DN and SXTs led to the suggestion that their 
outbursts have the same origin (van Paradijs & Verbunt 
1984; Cannizzo, Ghosh & Wheeler 1982; Huang & Wheeler 
1989; Mineshige & Wheeler 1989). 

Several important features of DN remain unexplained 
by the standard version of the DIM, which assumes that 
the disc extends down to the white dwarf's surface, that 
the mass transfer rate from the secondary star is constant, 
and that effects of disc and secondary irradiation can be 
neglected (see e.g. Lasota & Hameury 1998; Lasota 1999; 
Warner 1998; Smak 1999, Hameury. Lasota & Warner 1999b 
for a discussion of these effects). This version cannot ac- 
count for the quiescent accretion rates onto the white dwarf 
(they exceed the model predictions by more than two or- 
ders of magnitude), and for the very long recurrence times 
of WZ Sge-type systems (Lasota 1996a, b). The observed de- 
lay between the rise in optical and in EUV (e.g. in SS Cyg; 
Mauche 1996) was also long considered as a problem, but 
Smak (1998) recently argued that it can be naturally repro- 
duced by the DIM if appropriate outer boundary conditions 
for the disc are used. 

Similar difficulties, but more pronounced, appear for 
SXTs, since the standard DIM fails to reproduce a num- 
ber of essential properties of these systems. In particular, an 
early attempt to explain the properties of BH SXT outbursts 
with the DIM (Mineshige & Wheeler 1989) was unable to 
reproduce the observed long recurrence times (comparable 
to those of WZ Sge systems) and the high accretion rates re- 
quired in outburst. In addition, the weak X-ray flux detected 
from quiescent BH SXTs, which was mistakenly thought by 
some to be an argument in favour of the DIM versus the 
competing mass transfer instability model, is many orders 
of magnitude larger than any flux predicted by the DIM 
(Lasota 1996e). 

In the DN case, the problems can been solved, at least 
in part, if a "hole" is present in the inner regions of the 
disc. The "hole" can either be caused b y evaporation of the 
disc ( Meyer fc Meyer-Hofmeister 1994 ), or due to the ex- 
ten ded magnetospheric cavity of the W D (Livio & Pringle 
1992; Lasota, Kuulkers & Charles 1999). A model in which 
a "hole" is included can account for both the EUV delay 
and the quiescent X-ray luminosity (if the hole is filled with 
a tenuous. X-ray emitting gas; Meyer & Meyer-Hofmeister 
1994; Hameury, Lasota & Dubus 1999a). 

The case of WZ Sge requires additional modifications 
of the standard DIM. The very long recurrence time of this 
system, and the large amount of mass accreted during the 
outburst, can be reproduced either by using a value of the 
viscosity parameter acoid several orders of magnitude lower 
than in other DN (Smak 1993; Osaki 1995), or by making the 
"hole" sufficiently large for the remaining disc to be glob- 
ally stable. In the latter case, the outbursts would have to 
be triggered by enhanced mass transfer (Lasota, Hameury & 



Hure 1995; Hameury, Lasota & Hure 1997b). A marginally 
s table , truncated disc, as proposed by Warner, Livio & Tout 
( 1996| ) , i s equivalent to a globa lly stable disc from this point 



of view ( Hameury et al. 1997b). If the truncated disc is glob- 
ally stable, and QcoM ~ 0.02, however, there is not enough 
mass in the disc to account for the outburst's total energy, so 
that enhanced mass transfer is also required during the out- 
burst. Hameury et al. (1997b) find that if this enhanced mass 



transfer is due to irradiation of the secondary, the model re- 
produces the shape of the observed lightcurve. 

The presence of "holes" in the accretion discs of quies- 
cent SXTs seems unavoidable if their outbursts are due to 
the thermal- viscous instability (Lasota 1996a,b). In the BH 
case, such holes can only result from the evaporation of the 
accretion disc since black holes are not magnetized. Narayan, 
McClintock & Yi (1996) proposed a model in which, as a re- 
sult of evaporation, the inner parts of the accretion flow in 
quiescent BH SXTs form an advection-dominated accretion 
flow (ADAF; see also Lasota, Narayan & Yi 1996; Narayan, 
Barret & McClintock 1997). The model explains the ob- 
served spectra and solves the problem of incompatibility 
between the accretion rates predicted by the DIM and the 
observed flux. 

Esin, McClintock & Narayan (1997; see also Narayan 
1996; Esin et al. 1998) extended the idea of two-component 
accretion flows a step further. They showed that the various 
spectral states of black hole X-ray binaries can be modeled 
as a sequence of physical states with varying accretion rates 
and sizes of the two components of the accretion flow. Ac- 
cording to this model, in outburst, the disc extends down 
to the black hole, while the inner disc is gradually replaced 
by an ADAF as a system goes to lower luminosity levels. 
A gradual evaporation of the inner disc during the decline 
is also qualitatively consistent with the accretion geometry 
inferred by Zycki, Done & Smith (1999) from the X-ray re- 
processing properties of several BH X-ray transients. 

Hameury et al. (1997a) used the two component (inner 
ADAF + outer disc) model to describe the quiescent state 
and the rise to outburst of the BH SXT GRO J1655-40. 
The model predictions agree with existing constraints on 
the spectrum of the system in quiescence, the observed 6 
day delay between the rise to outburst of the optical and X- 
ray fluxes, and the rise-times of these two fluxes. Recently, 
Esin, Lasota & Hynes (1999) showed that the whole, unusual 
outburst can be described in the framework of the DIM thus 
validating Hameury et al. (1997a) results. 

In this paper, we consider the role of disc evaporation 
on the predictions of the DIM for SXT outbursts. We shall 
refer to a DIM that includes evaporation as a Truncated- 
Disc Instability Model (hereafter TDIM). Cannizzo (1998) 
already constructed TDIMs of the BH SXT A0620-00 with a 
weak evaporation and the prescription a = 50{H / R)^''^ for 
the viscosity parameter. Cannizzo's simulations show that 
models without evaporation lead to unreasonably short re- 
currence times, and that evaporation acts to increase the 
predicted recurrence times to values of several tens of years, 
as required. This is expected, since evaporation removes the 
inner, most unstable parts of the accretion disc. However, 
the outburst rise time predicted by Cannizzo's model is an 
order of magnitude longer than the observed one. This is due 
to the use of the prescription a = 50{H / R)^^^ according to 
Lasota & Hameury (1998). 
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There are two additional reasons to go beyond Can- 
nizzo's calculations. First, his calculations do not allow the 
disc outer radius to vary with time. Smak (1984, 1998; see 
also Hameury et al. 1998) showed that including this effect is 
crucial for reliable DIM predictions of DN outbursts, and it 
is likely to be important for SXT models as well. Second, in 
the models of Cannizzo (1998), the effect of including evap- 
oration and using a very small value of a in the quiescent 
disc are not clearly separated. In particular, in his work, it is 
unclear if evaporation alone can be held responsible for the 
long recurrence times, if values of a similar to those inferred 
from applications of the DIM to DN (typically ahot ~ 10~^ 
and acoid ~ 0.02) arc used for the discs of SXTs. In our view, 
this is an important question because the evidence for disc 
evaporation, even indirect, is stronger than for any specific 
a— prescription. 

In § 2, we recall some important observational proper- 
ties of BH and NS SXTs, later used to constrain our models. 
In § 3, we first describe the numerical truncated disc insta- 
bility model used. We then show that, in the limit where 
illumination can be neglected, TDIMs with values of a con- 
sidered as standard for DN discs are unable to reproduce 
the long recurrence times of SXTs, no matter how strong 
the evaporation is. In § 4, we present other plausible models 
of SXTs. We find that a TDIM with a small value of acow 
(~ 5 X 10~^) and strong evaporation does reproduce the 
long recurrence times and the high maximum values of the 
accretion rate in the discs of SXTs. We also present models 
in which the disc of SXTs are globally stable and the out- 
bursts are triggered by a slow variation of the mass transfer 
rate of the secondary. In § 5, we discuss important conse- 
quences, limitations and possible extensions of this work. In 
§ 6, we summarize our main results. 



2 MAIN OBSERVATIONAL 

CHARACTERISTICS OF SXTS 

BH and NS X-ray Transients show a rather complex be- 
haviour (see, e.g., Tanaka & Shibazaki 1996 and CSL for re- 
views). In this section, we identify several important proper- 
ties of these systems which may be considered as character- 
istic of their class and therefore must be explained by a suc- 
cessful model of SXT outbursts. We will also often compare 
SXTs with DN, since the same underlying thermal-viscous 
instability is supposed to be responsible for the outbursts in 
the two classes of systems. 

As a guide, we use the observed properties of SS Cyg, 
Aql X-1 and A0620-00 as indicative of the general proper- 
ties of U Gem-type DN, NS SXTs and BH SXTs respec- 
tively, keeping in mind that significant differences exist even 
between two outbursts in a same system (e.g. at least four 
types of outbursts in SS Cyg; Warner 1995). 



2.1 Timing properties 

The timescales of evolution of an unsteady disc are simple 
and robust predictions of the DIM, which provide powerful 
tests of the model when compared to the observations. This 
is because the time evolution of the disc is due to the exis- 
tence of a well defined limit cycle in the DIM, and because 



the observed variability is 'ready-to-interpret' (free of noise, 
systematic errors, etc.). 

Most of the properties of BH and NS SXTs mentioned 
in this section are taken from the compilation of CSL. We 
identify several important timescales for X-ray transients: 
the recurrence time Uec between two successive outbursts, 
the e-folding rise timescale Trise and the total rise time trise 
of an outburst (which can be wavelength-dependent), the 
e-folding decay timescale Tdecay (which is wavelength depen- 
dent), and the total duration of the outburst idm . In our 
view, a model which fails to reproduce all these character- 
istic times is not fully satisfactory. None of the previous 
models of SXTs passed this test (e.g. Mineshige & Wheeler 
1989; Carmizzo 1998). 

We focus our discussion on a sample of SXTs with 
clearly identified primaries (see Garcia et al. 1998): the BH 
systems have firm dynamical lower limits on the mass of the 
primary, and the NS SXTs showed type I X-ray bursts. We 
also focus on systems which have exhibited standard Fast 
Rise Exponential Decay (FRED) lightcurves. This excludes 
the long period systems GRO J1655-40 and V 404 Cyg 
and reduces the sample to four BH SXTs (A0620-00, GS 
2000+25, GRS 1124-683=Nova Mus 91 and GRO J0422+32 
- "best examples of FREDs" according to CSL) and two 
NS SXTs (Cen X-4 and Aql X-1; here we chose only their 
outbursts of the FRED type : it is not clear how the stan- 
dard DIM could account for the other types of lightcurves 
observed in SXTs). 

2.1.1 Recurrence times 

Aql X-1 experiences an outburst every year approximately, 
while the recurrence time of Cen X-4 is ~ 10 years. A0620- 
00 experienced two outbursts separated by 58 years (note 
that no X-ray data are available for the presumed first out- 
burst) . Most BH SXTs did not experience a second outburst 
since their discovery by X-ray satellites, so that there is good 
evidence for the recurrence times of BH SXTs generally ex- 
ceeding 20 — 30 years. As mentioned above, SXT recurrence 
times are similar to those of WZ Sge-type DN. For compar- 
ison, the recurrence time of SS Cyg is tree ^ 50 days (e.g. 
Cannizzo & Mattel 1992) and the recurrence time of other 
U Gem type DN is generally of the order of one month (e.g. 
Szkody & Mattel 1984). 

2.1.2 Rise times 

The characteristic X-ray rise timescales (rrise) of FRED- type 

SXTs vary from 0.3 to 2 days (for the period of fastest flux 
increase - CSL). The total rise time trise from quiescence to 
outburst peak is difficult to estimate because it is instrument 
dependent, and because the coverage of the rising phase is, 
for obvious reasons, rather poor. For SXTs, trise is ~ 5 to 
10 days in X-rays, or perhaps longer. For comparison, the 
total optical rise time in SS Cyg is triso ~ 3 days (for a 'type 
B' outburst - see below). The optical rise time of several U 
Gem- type DN in the sample of Szkody & Mattel (1984) is 
trise ~ 2 days. 

In several DN, there is evidence for at least two types 
of outbursts. 'Type A' outbursts correspond to FREDs and 
are usually interpreted as 'outside-in' outbursts, i.e. out- 
bursts triggered in the outer regions of the disc (Smak 1984). 
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This interpretation seems, in a few cases, consistent with the 
observed delay between the rise to outburst of optical and 
EUV (or UV) fluxes (Warner 1995). 'Type B' outbursts have 
a more symmetric shape with a slower rise phase, and are 
usually interpreted as 'inside-out' outbursts (Smak 1984; see 
also Menou et al. 1999a). 

In the case of SXTs, there is no observational evidence 
for outside-in outbursts. The delay between the rise of the 
optical and X-ray fluxes (observed only during a non-FRED 
outburst of the BIf SXT GRO J1655-40) is not an indication 
of an outside-in outburst if the disc is truncated. Indeed, 
Hameury et al. (1997a) obtain correct rise times and the X- 
ray/optical delay for an insidc-out outburst. In this case the 
'inside' corresponds to the inner edge of the truncated disc, 
far from the central object. 



2.1.3 Decay times 

The X-ray decay timcscalcs race of SXTs with FRED light- 
curves are typically 25 - 40 days (CSL). The optical decay 
timescales of BH SXTs are ~ 50 to 200 days, while the 
optical decay timescales of NS SXTs are ~ 10 to 30 days (but 
the available information is poor). We note that long decay 
timescales are also observed in WZ Sge-type DN (Kuulkers, 
Howell & van Paradijs 1996; Kuulkers 1998). 



2.1.4 Duration times 

The reported duration times for the FRED-type outbursts 
of BH SXTs range from 170 to 260 days in X-rays (and 
even longer in optical), while tdui ~ 70 days in X-rays for 
FRED-type outbursts of NS SXTs. The typical duration 
time of U Gem-type DN in the sample of Szkody & Mattei 
(1984) is tdur ~ 10 days (optical). 

One has to remember, however, that the reported du- 
rations for SXTs depend on the detection threshold of the 
instrument and that the coverage of some outbursts is in- 
complete; sometimes the reported duration is only a lower 
limit on the true duration. In some cases, the decay timescale 
Tdec may be more appropriate than the duration time tdec for 
comparison of the model predictions with the observations. 



2.2 Other important observational constraints 

An estimate of the maximum luminosity reached in outburst 
by a SXT depends on the distance to the system and, when 
compared to the Eddington luminosity, on the mass of the 
primary. According to Garcia et al. (1998), the NS SXTs 
considered here approximately reached the Eddington lu- 
minosity at outburst peak, as did some of the BH SXTs. 
According to CSL, the maximum luminosities reached are 
closer to 10% of the Eddington luminosity, but their val- 
ues could be underestimated. For example, CSL obtain for 
Nova Muscae 1991 a peak luminosity (0.4 - 10 keV) of 0.08 
in Eddington units, whereas according to Esin et al. (1997), 
who used the entire available spectral information, the peak 
luminosity in this system was almost exactly Eddington. In 
our view, this means that models of SXTs have to be able 
to reach the Eddington accretion rate. The Eddington ac- 
cretion rate is defined here as 



MEdd = 1.39 X 10'* (^-^^ g s-\ (1) 

the rate at which an accretion flow around a central object of 
mass Ml reaches the Eddington luminosity for a standard 
10% radiative efficiency. The Eddington luminosity, how- 
ever, seems to be the maximum luminosity for SXTs. Com- 
bined with outburst durations this implies that the mass 
accreted during outbursts by the central body is ~ a few 
lO'^^ g (roughly the same as accreted onto the white dwarf 
during WZ Sgc eruptions). This constraint must be satisfied 
by SXT models as well. 

Observations of broadened Ha emission lines provide 
an upper limit on the value of the disc inner radius i?in 
in quiescence for several BH SXTs (e.g. Orosz et al. 1994; 
Narayan et al. 1996; Menou, Narayan & Lasota 1999b). In 
most systems, the constraints indicate a value of _Rin less 
than 10"* to a few W^Rs, where Rs = 2.95x10^ {Mi/Mq) cm 
is the Schwarzschild radius of the black hole primary of mass 
Ml. The limits provide a useful constraint on the strength 
of evaporation to be included in the TDIMs of SXTs (see 
below). 

For some systems, other constraints on Rin exist at var- 
ious stages of the outburst from the modeling of X-ray re- 
flection properties (Zycki ct al. 1998), as well as spectral 
fits which involve Rin as a parameter (e.g. Sobczak et al. 
1999). However, these two methods are model dependent 
and somewhat controversial, because the relevant spectral 
features can be modeled differently (M. Nowak, private com- 
munication). We choose to ignore them here. 



3 "STANDARD" TRUNCATED DISC 
INSTABILITY MODELS OF SXTS 

In this section, we investigate a "standard" version of the 
TDIM for SXTs. As mentioned above, by "standard", we 
mean a model with values of the viscosity parameter a that 
are usually considered eis standard for DN discs: ahot ~ 0.1 
and Qfcoid ~ 0.02. Non-standard models, in particular models 
using significantly smaller values of Ocoid, are considered in 
§ 4. 

3.1 Numerical specifications 

3.1.1 General 

We use the numerical code described in Hameury et al. 
(1998) to simulate the evolution of discs around NSs and 
BHs. The equations of conservation of mass, angular mo- 
mentum and energy for a thin Keplerian disc are solved on 
an adaptive grid which resolves narrow structures in the disc 
(transition fronts in particular). A fully implicit numerical 
scheme is used for time evolution, to avoid any Courant-type 
limiting condition on the timestep of integration. 

A grid of disc vertical structures, which give the local 
cooling rate of the disc as a function of its surface density, 
central temperature and the vertical gravity is precalculatcd 
before running the disc evolution. The effect of changing 
the mass of the central object (from a NS SXT to a BH 
SXT) enters the time-dependent calculations both through 
the radial conservation equations solved, and through the 
grid of vertical structures with different gravities. 
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The outer radius of the disc is allowed to vary with time 
under the action of the torque applied by the secondary star 
on the outermost regions of the disc. Including this proper 
boundary condition (as opposed to a rigid boundary condi- 
tion with a fixed value of the outer radius) has been found 
to drastically affect the predictions of the DIM for DN discs 
(see Hameury et al. 1998 for details). The effect of including 
this proper boundary condition for discs of SXTs is discussed 
in § 3.2.1. 



3.1.2 Formation of a hole in the disc 

The formation of a hole in the discs of SXTs is proba- 
bly determined by the physics of the evaporation of gas 
from the disc, possibly into an ADAF. This physics is 
poorly understood. Several mechanisms have been proposed 
in the literature to explain the transition (see, e.g., Meyer & 
Meyer-Hofmeister 1994; Liu, Meyer, Meyer-Hofmeister 1997; 
Narayan & Yi 1995; Honma 1996; Shaviv, Wickramasinghe 
& Wehrse 1999), but none of these models has reached a 
sufficient level of detail to allow robust quantitative predic- 
tions. 

Disc evaporation has already been included in the cal- 
culations of Hameury et al. (1997a) and Hameury et al. 
(1999a). In these studies, a 'reasonable' formula was used 
(close to the one proposed by Meyer & Meyer-Hofmeister 
1994) and evaporation was included in the numerical mod- 
els as a sink of matter in the mass conservation equation 
(see also the models with evaporation of Cannizzo 1998). 

Here, we chose a slightly diff'erent approach. We assume 
an ad hoc law for the rate of evaporation as a function of 
radius, as in Hameury et al. (1997a; 1999a), but in addition 
we assume that evaporation acts only in a very narrow region 
at the inner edge of the disc. This is a good approximation 
if the rate of evaporation of the disc material is a steep 
function of radius. 

To be more specific, the inner radius of the disc, Rin, 
varies according to: 



Afdisc(-Rin) — Afcvap(-Rin) 



(2) 



where Mdisc(ftn) is the time-dependent accretion rate at the 
inner edge -Rin of the disc, and Mevap{R) is the evaporation 
law with an arbitrary functional form. This evaporation law 
is chosen so that, in outburst, Mcvap(-R) is much less than 
the accretion rate in the disc at any radius, and the disc 
is fully extended. The radius of truncation of the disc in 
quiescence depends on the functional form of M^vapiR) and 
the profile of Mdisc in the disc. We find that this numerical 
recipe is much more stable than that used in Hameury et al. 
(1997a), while it gives very similar results for the evolution 
of i?in and the outburst cycles in general. 

Following Hameury et al. (1997a), we could have 
assumed that the disc is evaporated into an ADAF 
which accretes everywhere at the maximum rate allowed, 
M^Jl^^ (R), and take for the evaporation law MevapiR) = 
M^Ii^^iR)- Recently, however, the estimates of M^'^^'^{R) 
have been revised to higher values (Esin et al. 1997; Menou 
et al. 1999b). Using these recent estimates for the evapora- 
tion law (i.e. a stronger evaporation than, e.g., in Hameury 
et al. 1997a) would imply that truncated discs are glob- 
ally stable (and stationary) for all mass transfer rates below 
~ 10^^ g s~^ (this is demonstrated by the Fig. 5 of Menou 



et al. 1999b), whereas typical SXT mass-transfer rates are 
much lower: ~ lO^'^ - 10^^ g s"^ (van Paradijs 1996; CSL; 
Menou et al. 1999b). In addition, for mass transfer rates 
< lO^'^ g s~^, the stable discs would be truncated at radii 
> 3 X 10* Rs, which is inconsistent with the upper limits 
derived from Ha emission lines (§ 2.2). 

To circumvent these difficulties, and ensure that discs in 
our models experience outbursts even at lower mass transfer 
rates, Af, 

following evaporation law 



cvap{R) must be less than M^^g^^(i?). We use the 



Mevap(-R) 



0.08AfE 



(3) 



where 5 is a constant ensuring that a maximally truncated 
disc is still unstable. The power laws and scalings in Eq. (3) 
are such that the profile of AIcvap{R) is initially slowly de- 
creasing with radius, and then decreases much more steeply 
beyond a radius ~ 8007?s x £~^^^. In all the models with 
evaporation described below, we take £ — 20 (except in sec- 
tion 4.1, in which £ — 30). This corresponds to a situation 
where the evaporation rate during quiescence is significantly 
smaller than the mass transfer rate from the secondary, thus 
ensuring that the disc is not marginally unstable. 

This expression for the evaporation law is ad hoc but 
possesses several advantages. First, it is a continuously de- 
creasing function of radius, which seems quite natural, and 
is appropriate for the numerical implementation of evap- 
oration described above. Second, the maximum value of 
Mcva.p{R), as defined in Eq. (^ when R ~ Rs, is consis- 
tent with the high value required in the ADAF models of 
Esin et al. (1997) to fit the high luminosity states of the 
BH SXT Nova Muscae 91. Third, it leads to disc truncation 
radii in quiescence which are smaller than the upper lim- 
its inferred from the broadening of Ha emission lines. The 
resulting discs are unsteady and experience limit cycles for 
reasonable mass transfer rates (~ lO'^^ g s~^), as required 
for the DIM to be the explanation of SXT outbursts. 

Naturally, one could worry that the predictions of the 
time-dependent models depend on the evaporation law used. 
This is true about the strength of evaporation, but we find 
that our conclusions do not crucially depend on the pre- 
cise shape of the evaporation law. We discuss the effects of 
changing the evaporation law on the model predictions in 
§5. 

Note that in the models presented here, the presence of 
a tenuous region in the inner accretion flow, presumably an 
ADAF, is only taken into account through disc evaporation 
(i.e. as a sink of mass): we neglect any other physical effect 
that this inner region of the accretion flow could have on the 
disc evolution. Disc irradiation by the central X-ray source 
during high luminosity states is also neglected, although it is 
a potentially important effect (see § 5 for a discussion). For 
simplicity, we use the same evaporation law for BH and NS 
SXTs, except that A^Edd and Rs are scaled with the actual 
primary mass Afi. 



3.2 Results 

We present in this section two SXT models, one with pa- 
rameters relevant for the BH SXT A0620-00, and the other 
for the NS SXT Aql X-1. In both models, ahot = 0.1, 
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Figure 1. Evolution of the disc mass Mdiso the V band magni- 
tude of the disc, the accretion rate at the disc inner edge Min, the 
disc outer radius -Rout and the disc inner radius Hin in a model 
with parameters relevant for the BH SXT A0620-00 (see text for 
details). The recurrence time is tree — 5 years in this model with 
aeold = 0.02. 
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Figure 2. Details on the evolution of the V band magnitude of 
the disc, the accretion rate at the disc inner edge and the 

disc inner radius ij;,; during an outburst of the model presented 
in Fig. 1. Note the multiple reflares during the decline, as seen in 
the modulations of the V band magnitude, Mjn and 



Qcoid = 0.02, the disc inner radius is limited during the 
evolution to the minimum value -Rmin = 5 x 10* cm (see §5 
for a discussion of this numerical limitation) and Eq. (pf), 
scaled with Mi, is used for the evaporation law. 

The model of A0620-00 has the following parameters: 
primary mass Mi = QMq, time-averaged value of the disc 
outer radius < Rout >~ 1.24 x 10^^ cm and mass trans- 
fer rate Mr = 3 x 10^^ g s~\ The model of Aql X-1 has 
the following parameters; primary mass Mi = 1.4Af0, time- 
averaged value of the disc outer radius < iiout >= 1-8 x 10^^ 
cm and mass transfer rate Mt = 2 x 10^*^ g s"\ The mam 
properties of these models are summarized in Figures |l| to ^ 
and are discussed in more detail below. 

For comparison, we also show in Fig. ^ and ^ results 
from a model with parameters relevant for A0620-00, but 
without evaporation. Results are similar, except for a signif- 
icantly shorter recurrence time, a larger number of reflares, 
and the existence of two type of outbursts. Successive small 
and large outbursts are also predicted in models of dwarf 
novae (see e.g. Cannizzo 1993a; Menou et al. 1999a). Dur- 
ing small outbursts, the outward propagating heating front 



does not reach the outer disc radius. In the DIM, the inner 
regions of quiescent accretion discs are more subject to in- 
stabilities than the outer ones, so that the truncation of the 
inner disc has a stabilizing effect. This explains the absence 
of 'small' outbursts and the reduced number of reflares in 
models that include evaporation (see below in § 3.2.6). 

The panel representing the total mass in Fig. ^ shows 
that the disc evolution is periodic and therefore relaxed (i.e. 
does not depend on the arbitrary properties of the initial 
disc used at the beginning of the simulation). This is true of 
every model discussed in this paper. Strictly periodic out- 
burst cycles are produced in SXT models, like for DN, while 
observed cycles show, in general, only some regularity. 

In the models presented here (as well as in all the other 
sections) , the outbursts are triggered in the inner regions of 
the truncated disc (inside-out type). Such outbursts never 
start exactly at the disc inner edge, so that there is always 
a heating front traveling some distance inwards. Its effect is 
seen as a 'spike' in the light-curve in Fig. 4, which lasts for 
less than one minute and is therefore not observable. In a 
truncated disc, these short-lived fronts are absent, because 
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Figure 3. Evolution of the total disc mass M^isci visual magni- 
tude of the disc and mass accretion rate onto the compact object 
Min in a model with the same parameters as the one shown in 
Fig. 1^, but without evaporation. The inner disc radius is fixed at 
5 X 10* cm. Here successive small and large outbursts are present. 
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Figure 4. Details of one of the outbursts of Fig. 
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Figure 5. Evolution of the total disc mass M^^^^, the V band 
magnitude of the disc, the accretion rate at the disc inner edge 
Mjn and the disc inner radius -Rin for two outbursts in a model 
with parameters relevant for the NS SXT Aql X-1. The recurrence 
time is troc — 1 year in this model with acold = 0.02. Here reflares 
are absent. 



the innermost regions where they would propagate no longer 
exist. Heating fronts only marginally reach the outer edge of 
the disc in the model with parameters relevant for A 0620-00. 
In the model with parameters relevant for Aql X-1, the disc 
is so extended that the heating fronts never go beyond one 
third of the disc outer radius. This would of course not be 
the case for shorter period systems such as GRO J0422-f 32. 

In the following discussion of the rise times and decay 
times predicted by the models, we will assume for simplicity 
that the X-ray emission from the system is proportional to 
the accretion rate Afin at the inner edge of the disc. The 
real situation is likely to be more complicated, in particu- 
lar if some of the X-rays observed are emitted by an ADAF 
whose luminosity is not simply proportional to M (in that 
case, the timescales could be even shorter because the radia- 
tive efficiency increases with M). Similarly, the V magnitude 
shown are not supposed to directly represent the V magni- 
tude observable from a system but only the intrinsic optical 
luminosity variations of the disc. In observed systems, the 
V band magnitude emission has also contributions from the 
secondary, from reprocessing of X-ray photons by the disc 
and possibly from an ADAF. We do not attempt to de- 
scribe all these processes here, and we concentrate on the 
basic properties of the models. 
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3.2.1 Recurrence times 

As for DN discs, it is important to include a correct bound- 
ary condition at the outer edge of the disc of SXTs. We 
find that the recurrence times predicted by models in which 
Rout is allowed to vary with time are typically a factor of a 
few (say 2-3) times shorter than those predicted if Rout is 
kept fixed. Indeed, when a heating front reaches -Rout , mass 
is spread over a more extended region if -Rout is allowed to 
vary. Consequently the critical surface density Emin below 
which a cooling front appears in the disc is reached by the 
disc earlier, less mass is accreted during the outburst, and 
overall it taJies less time for the disc to reach the state of 
critical surface density Emax at which the next outburst is 
triggered (Hameury et al. 1998) . Therefore, all existing mod- 
els of BH SXTs with fixed values of -Rout have systematically 
overestimated the recurrence times t^ac by a factor of a few 
when the heating front was able to reach the disc's outer 
edge. In the present paper, we only consider models with 
varying -Rout • 

The recurrence time in the model of A0620-00 is tree ~ 5 
years, i.e. shorter than the known recurrence time by more 
than one order of magnitude (§ 2.1.1). The recurrence time 
in the model of Aql X-1 is tree ~ 1 year, in very good agree- 
ment with the know recurrence time. 

3.2.2 quiescence 

During quiescence, all of the accreted mass is processed 
through the ADAF and is therefore evaporated, so that 
the mass accretion rate shown in figs. 1 and 5 is the ac- 
tual evaporation rate except when the disc radius reaches 
its minimum value during outbursts. As can be seen, the 
evaporation rate remains significantly smaller than the mass 
transfer rate from the secondary, varying between 0.15 and 
0.60 Mt in the case of the A0620-00 model, and between 
0.10 and 0.15 Mt in the Aql X-1 case. These are close to 
the maximum evaporation rate for which the disc is unstable 
but they do not imply fine tuning. 

3.2.3 Rise times 

In the model of A0620-00, once an outburst is triggered, it 
takes the disc ~ 9 days to reach the compact object (the time 
required to go from 5 x 10* cm to ~ 10^ cm, not included 
in our simulations, is small because of the short timescales 
in the inner disc, so that we neglect it). It also takes ~ 14 
days for the accretion rate at the disc inner edge to reach 
its maximum value. Overall, the total rise time in X-rays is 
therefore trise ~ 5 days, which is in reasonable agreement 
with the available information on FRED-type light-curves. 
The total rise time in optical is much shorter because the 
outbursts are triggered in a truncated disc, i.e. effectively 
in the outer regions of the accretion fiow. Our results are 
similar to those of Hameury et al. (1997a) who addressed 
in more details the issue of time delays between the rise in 
optical and X-ray. 

The rise timescale for the period of fastest flux increase 
is Triso — 4 days in X-rays. This is somewhat longer than 
the observed rise timescales of 0.3-2 days (§ 2.1.2). For com- 
pleteness, we also give the rise timescale in optical, Tiisc — 1 
hour, for which there is very little available observational 



data (GRO J1655-40 being the only system in which the 
rise in optical was observed; see Orosz et al. 1997). This 
timescale is approximately the thermal time at the disc in- 
ner edge, since this is where most of the optical light comes 
from at the beginning of an outburst. It is not easily compa- 
rable to the observations, because of light dilution by both 
the secondary and the ADAF. 

In the model of Aql X-1 shown in Fig. |^, it takes the 
disc ~ 3 days to reach the compact object once an outburst 
is triggered. The fact that this time is shorter than in the 
model of A0620-00 can be partly attributed to a less efficient 
evaporation in this model (because -MEdd and Rs, in Eq. j^, 
scale with the mass Mi of the primary). It also takes ~ 12 
days for the accretion rate at the disc inner edge to reach 
its maximum value. Overall, the total rise time is therefore 
trisc — 9 days, which again is in reasonable agreement with 
the available data on FRED-type light-curves. 

The total rise time in optical is not so short: triso — 9 
days. This is because the disc is less strongly truncated 
than in the model of A0620-00, which results in much more 
"inside-out/symmetric" outburst shapes (e.g. Smak 1984). 
In addition, because the disc is large, the inner parts do 
not contribute as much as in the case of A0620-00. The rise 
timescale for the period of fastest flux increase is Trise — 1 
days in X-rays, which is consistent with available observa- 
tional data. It is Trisc — 1.7 hours in optical. 



3.2.^ Decay and duration times 

In the model of Aql X-1, the total duration time of outbursts 
is tdur — 25 days. This is short compared to the duration ^ 
70 days reported in the literature for FRED-type outbursts. 
Again, we note that the outbursts produced in this model are 
not really of the FRED-type, but have a rather symmetric 
"triangle" shape. 

The decay timescale early in the decline is rdec — 3.5 
days in X-rays and Tdcc — 4 days in optical. The timescales 
become smaller later on because the light-curves steepen. 
These values are clearly shorter than the observed values of 
25-40 days in X-rays and 10-30 days in optical. 

In the model of A0620-00, the total duration time of an 
outburst is tdur — 55 days. This is quite short compared to 
the typical duration times for BH SXTs with FRED-type 
light-curves, ~ 170 — 260 days in X-rays and even more in 
optical. Note that fdur would be shorter if refiares did not 
occur (see below for a discussion of refiares). 

Estimating the decay timescales is complicated by the 
presence of refiares in this model. For completeness, we es- 
timate two types of rdec- The overall Tdoc corresponds to 
the decay timescale for the global light-curve, ignoring the 
"fiuctuations" induced by the refiares. It is rdec — 20 days 
in X-rays and rdec — 22 days in optical. The actual decay 
timescales, in between two successive refiares, are shorter: 
rdec — 9 days in X-rays and rdec — 11 days in optical. We 
find it difficult to compare these timescales with observed 
values, since observed light-curves do not resemble the light- 
curves shown in Fig. ^ (see § 5 for a discussion) 

One should keep in mind that we did not take into ac- 
count outer disc X-ray irradiation which will increase signif- 
icantly both the duration and the decay timescale (Dubus 
1999a,b). 
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3.2.5 Outburst amplitude 

In the model of A0620-00, the maximum accretion rate 
reached at the inner edge of the disc is substantially sub- 
Eddington (~ O.OlSMEdd; Fig. |). This is far too smaU. In 
the model of Aql X-1, on the other hand, the Eddington 
limit is reached without difficulty (Fig. 

The total amount of mass accreted during an outburst 
is also quite small in the model of A0620-00. This is not 
surprising because the recurrence time is more than 10 times 
shorter than observed, and not much mass is accumulated 
in quiescence. On the other hand, the total accreted mass 
is comparable to what is deduced from the observations in 
the case of Aql X-1. Again, this is not a surprise since the 
recurrence time is about right. 



3.2.6 Reflares 

A striking feature of the models described here is the possi- 
ble occurrence of reflares during the decline of an outburst 
(Fig. [^). These reflares correspond to multiple reflections 
of cooling and heating fronts during their propagation in the 
disc. For the specified set of parameters, we find that they 
are present in models with a BH primary, but absent when 
the primary is a NS. The number of reflares is reduced when 
the disc is truncated because increasing the inner disc radius 
gradually removes the regions of the disc where reflections 
occur. 

Reflares are in fact a natural outcome of the DIM. They 
have been observed in many simulations over the years, in 
particular in models of DN discs in which a same value 
Qhot ~ Qcoid was used (e.g. Smak 1984). In that case, a 
cooling front develops, propagates over some fraction of the 
disc radial extent, is reflected into a heating front which goes 
back all the way up to the disc outer edge. This is repeated 
indefinitely and leads to small amplitude variations of the 
disc luminosity. This failure of the model has been overcome 
by using a higher value of a for the hot ionized disc than 
for the cold neutral disc, which results in large amplitude 
outbursts as those observed. 

The fundamental reason why a reflare occurs is because 
the surface density just behind a cooling front reaches the 
critical value Emax sometime during its propagation. This 



happens typically when the front is located at Rh, 



< 

lO^'^ cm in our simulations. An outward propagating heating 
front then appears which destroys the inner cooling front as 
mass starts being accreted efficiently in the cooling region. 
Later on, this heating front is reflected when the surface 
density just behind the front reaches the critical value Emin. 
A cooling front develops which shuts off the outward trans- 
port of angular momentum at the origin of the heating front 
propagation. A sequence of heating and cooling front reflec- 
tions corresponds to the multiple reflares observed in the 
simulations. 

What was probably not previously realized is that the 
occurrence of reflares depends not only on the values of a 
but also on the mass of the central object, as shown by our 
simulations. We note that reflares were also found in the 
simulations of the decline phase of BH SXT outbursts of 
Cannizzo, Chen & Livio (1995). Cannizzo et al. interpreted 
the reflares as numerical artifacts. This is not the case in our 



models: reflares are a generic property of the DIM and the 
TDIM for some range of masses and viscosity parameters. 

One can rather easily see how the presence or absence 
of reflares depends on the mass of the accreting object. As 
noticed by Menou et al. (1999a) the surface density at the 
cooling front is 



E(i?t, 



(Ri, 



(4) 



where Eniin(i?) is the 'minimum' critical surface density in 
the disc. Vishniac (1997) emphasized that the cold, outer 
regions of the disc behind the cooling front are essen- 
tially frozen during the cooling front propagation. Therefore, 
E(_R) behind the front is constant with time and equal to 
-/^(Mi, a)Emin, where K depends on the viscosity parame- 
ter Q and the mass of the accreting object Mi . According to 
the disc models of Menou et al. (1999a), for a central mass 
Ml = 1.2 Mq, K 4, while for Mi = 7 Mq, K 6 ~ 7. 
The ratio of the 'maximum' to 'minimum' surface densities 
is (Hameury et al. 1998) 

Emax (Ocold = 0.02) 



0.1) 



6.4, 



(5) 



(Oh o 

SO that for Mi — 1.2 Mq the post-cooling-front surface den- 
sity i^Emin < Emax, and uo reflarcs are expected. For AIi — 
7 Mq , the post-cooling-front surface density ifEmin 5^ Emax 
and it is not surprising that reflares are present. 

To explain the dependence of K on the mass of the ac- 
creting object, one can use the cooling front model of Vish- 
niac & Wheeler (1996). Although some of the predictions of 
this model were found to be somewhat inaccurate by Menou 
et al. (1999a), we expect the scalings derived here to be valid. 
In this model the front speed can be written as: 

"front °^ ("^) ■ (6) 

From mass conservation the post-cooling front density can 
be written as (see e.g. Vishniac 1997) 

Vr (front) 



"front 



(7) 



where (front) is the speed of matter at the cooling front. 
Using Vishniac & Wheeler (1996) model, one obtains 

Epf 



cx Mr 



(8) 



This scaling, K oc Mi , is in good agreement with the 
numerical results of Menou et al. (1999a). 

3. 2. 7 Consequences 

Our models confirm the results of Cannizzo (1998) showing 
that a stronger evaporation increases tree because it removes 
the inner, most unstable parts of the disc. The recurrence 
times are also increased if the mass transfer rate Mt in 
the model is reduced, since more time is then needed to 
accumulate mass in the disc up to the critical surface density 
Emax at which the next outburst is triggered. An exploration 
of the parameter space of the models shows that it is not 
possible, however, to reproduce recurrence times of ten to 
tens of years in TDIMs with 'standard' values of a. This is 
because to obtain these long recurrence times, the increase in 
the strength of evaporation or the reduction of Mr required 
lead to globally stable discs. 
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Figure 6. Recurrence times predicted by models of NS and BH 
SXTs (assuming parameters relevant to Aql X-1 and A0620-00, 
respectively) for various values of the disc fixed inner radius 
Rin (ahot = 0.1, Ocold = 0.02). The dotted, dashed and solid 
lines correspond to mass transfer rates Mt = 101*^,3 x and 
10^^ g s~^, respectively. For both NS and BH SXTs, truncating 
more the disc results in globally stable discs (tree — > oo) instead of 
tree > 10 — 50 years as expected. Models with arbitrary long t^cc 
can be constructed by fine tuning a marginally stable disc (e.g. 
the two upper points in the BH SXT panel), but these models do 
not have any physical meaning (see text for details). 



To demonstrate this, we computed a series of models 
of NS and BH SXTs with fixed values of i?in. The vari- 
ous values of iiin are equivalent to various strengths of disc 
evaporation. For example, a model with parameters relevant 
for A0620-00, with a fixed value of R\n = 6 x 10^ cm and 
Mt ~ 10^'' g s^^ has a recurrence time free — 3 years. This is 
slightly smaller than the recurrence time in the model with 
evaporation shown in Fig. |l|, in which Rin — 6 x 10^ cm 
before the onset of an outburst, but A4t is slightly smaller. 
This type of comparison shows that the main effect of evap- 
oration on the outburst cycles is to affect the radial extent 
of the disc in quiescence, which in turn determines tree- 

Figure ^ shows the recurrence times predicted for three 
mass transfer rates (M = 10^*^, 3 x 10^" and 10^^ g s"^) 
and several values of the fixed inner radius Rin- At a given 
Mt, the same trend is observed: the recurrence time initially 
increases rather slowly with iiin, but then the disc rapidly 
becomes globally stable at some critical value of Rin- For 
reasonable values of Mt, this critical value of iiin corre- 
sponds to recurrence times well below ten to tens of years. 
It is possible to adjust the inner disc radius so that the re- 
currence time will be arbitrary lorig (see e.g. the two upper 
points in the lower panel of Fig. |6|) but the resulting disc 
is marginally stable. Keeping it in this state for the entire 
quiescence period requires a mass-transfer rate Mt strictly 
constant, which is clearly unrealistic. 



This result on short recurrence times is confirmed by 
an analytical estimate of the recurrence times of SXTs. 

An outburst starts when the surface density reaches the 
critical value Emax. In the case of outside-in (type A) out- 
bursts, this happens when enough mass has accumulated at 
the disc outer edge. If the accumulation time at the outer 
edge is longer than the viscous time (as it is the case for low 
accretion rates), matter transferred from the secondary dif- 
fuses inwards and accumulates at shorter radii, giving rise 
to an inside-out (type B) outburst. Since, for a given set 
of parameters, the recurrence time of a type A outburst is 
shorter than the corresponding time for a type B outburst 
{tA < ts), we consider here only type B outbursts (in any 
case the only type obtained in our models) . 

Following Smak (1993), the recurrence time of a type B 
outburst is 



tB 



V dt ) 



(9) 



which is the characteristic growth-time of a surface density 
contrast. From the standard diffusion equation describing 
the evolution of S, one obtains the relation between ts and 
i/E, where v is the kinematic viscosity coefficient. The rela- 
tion between z^E and the effective temperature Tea is 

lyssT^s = g ^3 i^E, (10) 

where asB is the Stefan-Boltzmann constant, G is the grav- 
itational constant, R is the distance from the central object 
of mass Ml, and ^ ~ 3 — 5 is a constant factor that accounts 
for the non-stationarity of the quiescent disc (Idan et al . 
1999). Using the values of Emax calculated by Hameury et 
al. (1998) and assuming a constant effective temperature, 
which is a good approximation in quiescence (Smak 1984; 
Hameury et al. 1998), one gets 



vo.o2y uoook; 



yr- 



(11) 



This characteristic time ts is expected to be an upper limit 
to the recurrence time (e.g. Teff may be larger than 3000 
K). Again, this rough estimate is much less than tens of 
years. It also suggests that one way to obtain long enough 
recurrence times is to reduce the value of QcoM . Interestingly, 
we expect the maximum accretion rate reached by the disc 
during outburst to increase as well if Ocoid is reduced. This 
is because more mass will be present in the disc at the time 
an outburst is triggered (Emax oc a^'j'tf^; e.g. Hameury et 
al. 1998). Models of SXTs with smaller values of ftcoid are 
considered in the following section. 



4 "NON-STANDARD" TRUNCATED DISC 
MODELS OF SXTS 

In the following, we present two non-standard truncated disc 
models of SXTs. The first model is a non-standard TDIM 
because the value of the viscosity parameter cicoid used is 
significantly smaller than the standard value ~ 0.02. The 
second model explores the possibility that quiescent, trun- 
cated discs in SXTs are globally stable. This model is not a 
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Disk Instability Model, since the origin of an outburst has 
to be a variation of the mass transfer rate in the system. 

4.1 TDIMs with a smaller a in quiescence 

The model described in this section is similar to the model 
with parameters relevant for the BH SXT A0620-00 pre- 
sented in § 3, except for QcoW which is set to 5 x 10""^ here. 
We have also taken a slightly larger £ = 30. The predictions 
for the outburst cycles are shown in Fig. and a details on 
the evolution of important quantities during an outburst are 
shown in Fig. ^. 

As expected, the recurrence time predicted in this 
model is much longer {tree ~ 53 years) and the maximum ac- 
cretion rate reached in outburst is much larger (~ O.lMEdd) 
than in the model with acoid = 0.02. In our view, this is 
a success of the model that such long troc and high M in 
outburst can be reached with smaller values of acoid. 

There are other important differences caused by the re- 
duction of Qcoid- The disc mass is much larger in this model, 
as is the total mass accreted during an outburst. The varia- 
tions of -Rout are also larger. Finally, the disc optical emission 
in quiescence is quite reduced compared to what it is in the 
model with Ocoid = 0.02. 

The most obvious difference, however, is probably the 
disappearance of reflares in the decline phase of the out- 
bursts when ttcoid = 5 X 10~^ . This is in fact consistent with 
the explanation for the occurrence of reflares given in § 3.2.6. 
The ratio Emax/Smin is indeed increased from 7 to ~ 20 
when Qcoid is reduced from 0.02 to 5 x 10~^ (for the same 
Qhot), which makes it more difficult for the post-cooling- 
front surface density to reach Emax and trigger a reflare. 

The total duration of an outburst predicted in the 
model with Ocoid = 5 x 10~^ is tdur — 50 days (similar in op- 
tical and X-rays). This is the same as in the case Ocoid = 0.02 
and is still fairly short compared to the observed duration 
times (§ 2.1.4). 

Once an outburst is triggered, it takes the disc ~ 2 days 
to reach the compact object. The rise phase of the outburst 
(both in X-rays and in optical) can be clearly separated in 
an early phase of very rapid increase and a late phase of 
much slower increase (see Fig. |^ . The total rise time for the 
early phase in X-rays (up to Min ~ 2 x 10^ g s"'^; see Fig. ^ 
is trisc — 3 days, while the slower phase takes an additional 
~ 15 days to reach the outburst peak. The total rise time 
for the early phase in optical (up to V — 13; see Fig. ^ 
is triae — 1 day, while the slower phase takes an additional 
~ 10 days to reach the outburst peak. 

The rise timescales (we only consider the early phases 
of rapid increase here) are Trisc — 0.4 days in X-rays (short 
but consistent with observed values), and Triso — 0.1 days in 
optical. The decay timescales (early in the decline phase, be- 
fore it steepens) are rdoc — 10 days in X-rays and rdcc — 13 
days in optical. These values are too small, and not consis- 
tent with the lowest observed values. 

4.2 Models with globally stable truncated discs 

For completeness, we also constructed a model with param- 
eters relevant for the BH SXT A0620-00 in which the disc is 
truncated (same evaporation law as before), the mass trans- 
fer rate Afr is low enough for the disc to be globally stable 




50 100 150 

time (years) 

Figure 7. Evolution of the disc mass A/disc i the V band magni- 
tude of the disc, the accretion rate at the disc inner edge Min, the 
disc outer radius -Rout and the disc inner radius Hjn in a model 
with parameters relevant for the BH SXT A0620-00. The recur- 
rence time is tree — 50 years in this model with acold = 5 X 10~^. 

during a long 'quiescent' period, and an outburst is trig- 
gered by a slow variation of Mt- (Note that this is different 
from a burst of mass, or an enhanced mass transfer during 
the outburst.) More specifically, Mt was initially set to the 
value 10^"" g s~^ and was linearly increased over a period of 
several years until it reaches a value sufficient to trigger an 
outburst in the disc. The values of a chosen for this model 
are ahot =0.1 and QcoM ~ 0.02, i.e. similar to the 'standard' 
values of the models described in § 3.2. 

We find that in most respects, the predictions of this 
model for the outburst properties are similar to those of the 
TDIM of A0620-00 described in § 3.2 (which is why we do 
not show a figure for an outburst in this model; Fig. |^ is 
basically applicable): the maximum accretion rate reached 
is slightly smaller, and two reflares instead of three are pre- 
dicted (as well as a duration a bit shorter). 

Of course, there is no prediction for the recurrence time 
in this model, since one has to refer to some activity cycle 
in the secondary to obtain outburst cycles for the disc. 
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Figure 8. Details on the evolution of the V band magnitude of 
the disc, the accretion rate at the disc inner edge Af;,! and the 
disc inner radius iJin during an outburst of the model presented 
in Fig. 7 

. Note the absence of refiares. 



5 DISCUSSION 

We have constructed models of SXT outbursts that include 
evaporation of the disc inner regions during quiescence. The 
only numerical limitation in our calculations is the minimum 
radius -Rmin = 5 x 10** cm allowed for the disc inner radius 
-Rin during evolution. Our working assumption is that, in 
outburst, the disc extends down to the compact object with 
a nearly constant accretion rate in its inner regions. It seems 
a reasonable assumption given the results of existing calcu- 
lations on such extended discs (e.g. Cannizzo et al. 1995), 
and it allows us not to worry about modeling the innermost 
regions of the disc during the outbursts. 

One of the advantage of this simplification is that we 
do not have to handle the problematic radiation-pressure 
dominated regions of the disc, located close to the compact 
object in outburst. Our results concerning the recurrence 
times should not be affected by the numerical limitation 
-Rin > -Rmin siuce it doBS uot affect the propagation of the 
coohng wave and is irrelevant in quiescence (during which 

-Rin ^ ^^min)- 

There are, however, important physical effects that are 



not included in our simulations. Disc irradiation by the cen- 
tral X-ray source is not taken into account. Irradiation is 
likely to increase the duration times of the outbursts, as 
well as increase the recurrence times and the maximum ac- 
cretion rates reached in outburst, as suggested by King & 
Ritter (1998) and confirmed by numerical calculations of 
Dubus (1999a,b). Irradiation is also likely to increase the 
decay timescales, both in X-rays and optical. If the optical 
emission from the system is dominated by reprocessing of 
X-rays by the disc during outburst, the predicted optical 
light-curves could be markedly different from those found 
here in models without irradiation. Nevertheless, it is possi- 
ble that irradiation is unimportant in some SXTs (e.g. Shah- 
baz, Charles & King 1998). Our results should be applicable 
to these systems. 

Another effect that has been neglected in our simula- 
tions is the fact that the disc could reach the 3:1 resonance 
during its expansion in outburst, at least in short orbital pe- 
riod systems, where heating fronts can reach the outer edge 
of the disc. Tidal effects could then significantly affect the 
outburst properties (Osaki 1996). Similarly, enhanced mass 
transfer due to the irradiation of the secondary might play 
an important role in SXTs, as it has been argued for WZ 
Sge-type DN (Smak 1993; Hameury et al. 1997b; Augusteijn, 
Kuulkers & Shaham 1993). 

One of the advantages of our calculations is that they 
are free of all these complications, and isolate the effects of 
including the evaporation of the disc and of reducing the 
value of Ocoid- We acknowledge that none of the models 
presented in this paper matches the requirements that we 
defined in § 2 for a fully successful model of SXTs. This 
suggests that at least one ingredient of the model is miss- 
ing (the most obvious candidate being irradiation) . In future 
applications of the DIM to SXTs, additional physical effects 
will have to be included if ones wishes to reach a better 
agreement between the model predictions and the observa- 
tions. 

One of the predictions of the DIM that was clearly iden- 
tified in our models is the possible occurrence of refiares dur- 
ing the decline from outburst. Some of the BH SXT light- 
curves show 'refiares' during decline from outburst but these 
observed features are totally different from the refiares ap- 
pearing in our models. According to CSL, refiares (which 
they call 'secondary maxima') observed in X-ray transients 
are of three morphological types: 'glitches', 'bumps' and 
'mini-outbursts'. Glitches are upward infiections superposed 
on a smooth exponential decay. Neither the shape, nor the 
amplitude of refiares obtained in our models correspond to 
glitches. Bumps and mini-outbursts are more heterogeneous 
classes of features but it is difficult to find in the CSL com- 
pilation a light-curve morphology similar to the one seen 
in Fig. (Eh. The refiares are only present in models with 
a 3> 10 , i.e. in models with recurrence times too short 
to correspond to the observations. It is, therefore, a suc- 
cess of the model that the set of parameters ensuring both 
long recurrence times and absence of unobserved refiares are 
consistent. Some of the observed refiares may be due to disc 
irradiation (G. Dubus, private comm. 1999), which means 
that this effect has to be included before any detailed model 
of SXT with refiares can be constructed. 

The fact that, in our simulations, refiares occur orfiy 
when the mass of the accreting compact object is several 
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Mq (the boundary being ~ SM© for standard values of a) is 
interesting because it is apparently consistent with the fact 
that reflexes have been observed only in SXTs containing 
BHs, and not in SXTs containing NSs (e.g. CSL). However, 
since our roflarcs arc different from the observed ones and 
absent in the favoured model, this might be a coincidence. 
On the other hand, since a complete, self-consistent model of 
X-ray transients is still to be constructed it could be useful 
to keep in mind this property differentiating neutron star 
from black holes. 

By using the strongest evaporation compatible with ex- 
isting upper limits on the value of R\n in quiescent BH SXTs, 
we tested if evaporation alone can be responsible for the long 
recurrence times of SXTs, or if very small values of QcoM in 
their quiescent discs are also required. Our calculations in 
§ 3.2.7 show that a strong evaporation is not sufficient, and 
that either a viscosity parameter acoid at the level of ~ 10~^ 
is required to reach agreement with recurrence times of tens 
of years, or that disc outbursts must be drastically affected 
by X-ray illumination of the disc. 

The very low values of acoid required to model BH 
SXTs would suggest that more than one viscosity mecha- 
nism drives accretion in the disc. The viscosity in the hot 
state arises probably from a well developed MHD turbu- 
lence with a universal viscosity parameter Ohot ^0.1. The 
much lower value of a in quiescent discs is consistent with 
the disappearance of MHD turbulence in quiescent discs pro- 
posed by Gammie & Menou (1998) (see also Meyer & Mcyer- 
Hofmeister 1999). We note that if the relevant mechanisms 
in quiescence were non-local in nature, the use of a local 
a— prescription would be questionable and a major revision 
of the DIM might be required. 



6 CONCLUSION 

In this paper, we investigated the effect of disc evaporation 
on the predictions of the DIM for SXTs. As for DN discs, we 
show that the effect of allowing the disc outer radius Ro-at 
to vary with time is crucial to obtain robust predictions for 
the outburst cycles. 

We use the strongest evaporation still compatible with 
available upper limits on the value of the disc inner radius 
in quiescent BH SXTs to test if evaporation alone can be 
responsible for the long recurrence times of SXTs. 

We find that such a strong evaporation does increase 
the recurrence times predicted (in agreement with earlier 
calculations by Cannizzo 1998), but is not sufficient to re- 
produce the longest known recurrence times of SXTs when 
standard values of the viscosity parameter a are used in the 
disc (tthot ~ 0.1, acoid ~ 0.02). 

We show that models with strong evaporation and a 
significantly smaller value of acoid (- a few 10"^) do repro- 
duce the long recurrence times of tens of years of BH SXTs, 
as well as the high X-ray luminosities reached by these sys- 
tems in outburst. The value of acoid needed to reproduce 
the long recurrence times would be even smaller if evapora- 
tion wEis weaker than what we assumed. Another possibility 
is of course that X-ray illumination drastically alters the 
outburst cycle. We defer the examination of this effect to a 
future work. 

We argue that the requirement for a smaller value of 



Qcoid in at least some SXTs is consistent with the disap- 
pearance of MHD turbulence in quiescent discs (Gammie 
& Menou 1998) and suggests that another mechanism, per- 
haps non-local, could be responsible for accretion in the disc 
during this low luminosity phase. 

Finally, since the version of the DIM used in this pa- 
per fails to reproduce several important properties of SXTs 
(such as duration and shape of the light-curve) , future work 
must include in the model the observed missing ingredient: 
X-ray irradiation of the disc. 
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